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We have used normal metal-insulator-superconductor (NIS) tunnel junction pairs, known as SINIS struc- 
tures, for ultrasensitive thermometry at sub-Kelvin temperatures. With the help of these thermometers, we 
have developed an ac-technique to measure the electron-phonon (e-p) scattering rate directly, without any 
other material or geometry dependent parameters, based on overheating the electron gas. The technique is 
based on Joule heating the electrons in the frequency range DC- 10 MHz, and measuring the electron tem- 
perature in DC. Because of the nonlinearity of the electron-phonon coupling with respect to temperature, 
even the DC response will be affected, when the heating frequency reaches the natural cut-off determined 
by the e-p scattering rate. Results on thin Cu films show a T 4 behavior for the scattering rate, in agreement 
with indirect measurement of similar samples and numerical modeling of the non-linear response. 



1 Introduction 

Interaction between conduction electrons and thermal phonons is elementary for many processes and phe- 
nomena at low temperatures, and better measurements of the scattering rate need to be developed to char- 
acterise these processes accurately. Most earlier data was taken at DC, where sample parameters have to 
be taken into account before one can obtain the strength of the e-p interaction from a fitting parameter 
0] El ID . It would therefore be quite benefitial if one could measure the e-p scattering rate 1/t c _ p di- 
rectly without any fitting parameters and as a function of any external parameter. Some steps towards this 
direction have already been taken with a RF-based technique 1 5 1. 

Here we report a technique, which measures the e-p scattering rate directly without any fitting parame- 
ters and without the need of more complex RF circuitry. The rate is read directly from a cut-off frequency 
seen in the non-linear response of the DC electron temperature, measured with the help of NIS tunnel 
junctions at sub-Kelvin temperatures. 

In the disordered limit ql << 1, where q is wavevector of the dominant thermal phonons and / is the 
electron mean free path, theory predicts that the electron-phonon scattering rate from acoustic phonons is 
l/ r c-p = a ^e El El- where a is a sample parameter and T e the electron temperature. In contrast, in the 
pure limit temperature dependence 1/t c _ p = a'T e 3 is expected |8|. 
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2 Experimental techniques 



Our sample geometry is nearly identical to the one reported for DC measurements J4][9), schematic of the 
sample and the measurement circuit is depicted in Fig. ^a)- The sample used in this work has two Cu 
normal metal islands of length ~ 500/itn, width 300nm, and thickness 45nm, on oxidized Si substrate. 
The islands were separated by a distance of 2/im. The two aluminum wires going up in Fig. ^a) from 
each island were separated from Cu by a thin oxide layer, forming normal metal-insulator-superconductor 
(NIS) tunnel junctions in the overlap area. NIS junctions (SINIS structure) were used to measure the 
electron temperature of the Cu islands. Two aluminum wires going down from the lower Cu line were 
in metallic contact with the island, forming normal metal-superconductor (NS) junctions. These provided 
good electrical contact to the Cu, but the heat leak through them is negligible due to Andreev reflection, 
thus providing uniform heating power for the Cu wire. Heat leak through the tunnel junctions was also 
estimated to be at least an order of magnitude smaller than the cooling by electron-phonon coupling. 

SINIS thermometers were calibrated against a Ruthenium Oxide thermometer attached to the sample 
stage. The SINIS thermometers were current biased, and temperature of the cryostat (base temperature 
60 mK) was changed slowly to keep the sample in thermal equilibrium with the sample stage, while the 
voltages of the SINIS thermometers were measured. There was a small discrepancy with the measured cal- 
ibration curve and the corresponding BCS theory result at the lowest temperatures, but it can be explained 
by extra noise coming from the measurement circuitry to the sample l4l|9j. 
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Fig. 1 (a) The measurement setup used. Black horizontal lines describe copper, and grey vertical lines aluminum, (b) 
The measured DC (average) electron temperature vs. heating frequency. Each line is an average of about 20 frequency 
sweeps. 



The experiment then proceeds by heating the electrons of the lower wire via the NS junctions with an 
ac-voltage. The frequency of the heating signal was changed step by step, and the voltages of both SINIS 
structures were measured simultaneously. The upper Cu island in Fig. ^a) could be used to measure the 
phonon temperature |4, 9|, whereas the lower SINIS measured the heated electron temperature. We took 
many frequency sweeps with the same ac-voltage amplitude, and averaged the data to reduce the noise. The 
amplitude of the ac-voltage was then changed to obtain a different average electron temperature. A typical 
sweep consisted of 300-500 frequency points between 0. lHz-5MHz. The five representative sweeps shown 
in Fig. ^b) are an average of approximately 20 frequency sweeps each. Overheating is clearly seen in the 
high frequency range between lOkHz-lMHz for all the sweeps shown, which can be understood by the 
numerical modeling discussed below. 
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3 Numerical modeling 

To fully understand the effect of ac heating at a signal frequency uo on the electron temperature T e , we need 
to solve the full differential equation governing the heat flow: 



T e ^j- = —A(T™ - T p ") + P'[l - cos(2wt)], (1) 

where we have defined A = £'/7 and P' = P ac /(2'yil), where £' is the sample constant in the 
electron-phonon interaction power P = £'f2(T™ — T"), fl is the sample volume, 7 is the Sommerfeld 
constant and P ac is the amplitude of the ac power. Since £' also depends linearly on 7, we point out that 
the temperature relaxation rate does not depend on the heat capacity at all, and is only dependent on the 
electron-phonon scattering rate. 

We have integrated Eq. Q numerically for exponents n = 5,6, corresponding to 1/t c _ p = aT™, 
m = 3,4. From the T e (t) curve, the steady state average electron temperature has been determined as a 
function of u>, with varying ratio of ac heating P ac to DC heating power E'SIX^ 1 , and with realistic sample 
parameters. Fig. |2ja) shows the obtained results for four different power ratios in the case n = 6 (the 
impure limit). As can be seen, the average T e develops a clear step up at some cut-off frequency, and this 
frequency moves up as a function of T e . To be able to define this cut-off unequivocally, we analyse the 
data further by taking the logarithmic derivative with respect to cj, shown in Fig. EJb). The log-derivative 
develops a clear peak, whose position we will use as the definition of the cut-off frequency ui c . In addition, 
the log-derivative shows quite clearly the effect of the strength of the non-linearity: the larger the relative 
ac power (top curve largest), the more non-symmetric the peak is. 

If one plots ut c vs. the low-frequency limit of T e , one obtains the plot shown in Fig. |3ja), where both 
cases n = 5 and n = 6 are shown. It is quite clear that u c follows the right temperature dependence 
expected for a cut-off determined by the e-p scattering rate only, as long as the nonlinearity is not too 
severe (the lowest temperatures deviate a bit in this case). Therefore, we can state with confidence that the 
analysis scheme outlined above will be a direct measurement of the e-p scattering rate without any fitting 
parameters. 




Fig. 2 (a) Simulated average electron temperature, (b) Logarithmic derivative of the simulation data. 
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Fig. 3 (a) Maximum of the logarithmic derivative of the simulated data (squares) vs. T e . lj* = nAT?-' 2 . (b) 
Maximum of the logarithmic derivative obtained from the measured data (squares). Dotted line correspond to f c cx T 3 
and dashed to f c oc T 4 . Solid line is a best fit f c — aT n , with n = 3.7. 

4 Experimental results 

Figure|5Jb) shows the experimental data analyzed as outlined above, with the maxima of the log-derivatives 
of the measured data vs. the low-frequency limit of T e . The best two parameter fit f c — aT n gives n = 3.7. 
However, the scatter in the data is still a bit high, and to make a more qualitative analysis, we also compare 
the integer powers n — 3 and n = 4. As can be seen, the dashed line corresponding to power T 4 fits to 
data much better than the one corresponding to T 3 . This temperature dependence is in agreement with the 
DC measurements in the disordered regime |4][9). From the T 4 fit we get S' = 1.0 X 10 10 W/K 6 m 3 , using 
a literature value 7 = 0.69 m "^, 2 for copper, again in agreement with our DC data. One should keep in 
mind, however, that the numerical value of £' is strongly dependent on the sample disorder. In fact, we 
have seen in other, very similar samples that the value of £' can be almost an order of magnitude smaller. 

5 Conclusions 

We have presented a method to measure the electron-phonon scattering rate directly without the need to 
know any material or geometry dependent sample parameters. The measurements can be performed at DC, 
greatly simplifying the electronics, due to the non-linearity of the e-p interaction. Our initial experimental 
results agree with the DC data for Cu thin films in the disordered limit. 
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